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Abstract 

The three domatic number problem asks whether a given undirected graph can be 
partitioned into at least three dominating sets, i.e., sets whose closed neighborhood 
equals the vertex set of the graph. Since this problem is NP-complete, no 
polynomial-time algorithm is known for it. The naive deterministic algorithm for 
this problem runs in time 3", up to polynomial factors. In this paper, we design an 
exact deterministic algorithm for this problem running in time 2.9416™. Thus, our 
algorithm can handle problem instances of larger size than the naive algorithm in 
the same amount of time. We also present another deterministic and a randomized 
algorithm for this problem that both have an even better performance for graphs 
with small maximum degree. 

Key words: Exact algorithms, domatic number problem 

1 Introduction 

In this paper, we design a deterministic algorithm for the three domatic number 
problem, which is one of the standard NP-complete problems, see Garey and 
Johnson |GJ79|. This problem asks, given an undirected graph G, whether or not 
the vertex set of G can be partitioned into three dominating sets. A dominating set is 
a subset of the vertex set that "dominates" the graph in that its closed neighborhood 
covers the entire graph. Motivated by the tasks of distributing resources in a computer 
network and of locating facilities in a communication network, this problem and the 
related problem of finding a minimum dominating set in a given graph have been 
thoroughly studied, see, e.g., ICH77IFar84IBon85IKS94IHT98IFHK00IRR04l . 
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The exact (i.e., deterministic) algorithm designed in this paper runs in exponential 
time. However, its running time is better than that of the naive exact algorithm for 
this problem. That is, we improve the trivial 0(3") time bound to a time bound 
of (5(2.9416"), where the notation neglects polynomial factors as is common for 
exponential-time algorithms. The point of such an improvement is that a 0(c n ) 
algorithm, where c < 3 is a constant, can deal with larger instances than the trivial 
0(3") algorithm in the same amount of time before the exponential growth rate 
eventually hits and the running time becomes infeasible. For example, if c = V3 « 

1.732 then we have (5 = 0(3 n ), so one can deal with inputs twice as large as 

before. Doubling the size of inputs that can be handled by some algorithm can make 
quite a difference in practice. 

Exact exponential-time algorithms with improved running times have been 
design ed for var ious other important NP-complete problems. For example, Dantsin 
et al. |DGH + 02| pushed the trivial (5(2") bound for the three satisfiability problem 
down to 0(1.481"), which was further improved to (5(1.473™) by Brueggemann 
and Kern IBK04I . Schoning ISch02l . Hofmeister et al. [HSSW02| and Paturi et 
al. | PPSZ98 1 proposed even better randomized algorithms for the satisfiability problem. 
Combining their ideas, the currently best randomized algorithm for this problem is due 
to Iwama and Tamaki |IT03 1, who achieve a time bound of (5(1.324™). 

The currently best exact time bound of (5(1.211™) for the independent set problem 
is due to Robson HRob86l . Eppstein [Epp0 1a|Epp01b| achieved a (5(2.415") time 
bound for graph coloring and a (5(1.3289") for the special case of graph three 
colorability. Fomin, Kratsch, and Woeginger |FKW04| improved the trivial (5(2") 
bound for the dominating set problem to 0(1.93782"). Comprehensive surveys on 
this subject have been written by Woeginger | Woe03 1 and Schoning | Sch05 1 . 

In designing domatic number algorithms, it might be tempting to exploit known 
results (such as Eppstein's 0(1.3289") bound) for the graph three colorability problem, 
which resembles the three domatic number problem in that both are partitioning 
problems. However, as Cockayne and Hedetniemi ICH77 I point out, the theory of 
domination is dual to the theory of coloring in the following sense. Coloring is based 
on the hereditary property of independence. A graph property is hereditary if whenever 
some set of vertices has the property then so does every subset of it. In contrast, 
domination is an expanding property in that every superset of a dominating set also 
is a dominating set of the graph. Further, graph colorability is a minimum problem, 
whereas the domatic number problem is a maximum problem. Independence (and thus 
colorability) can be seen as a local property, since it suffices to check the immediate 
neighborhood of a set of vertices to determine whether or not it is independent. In 
contrast, dominance is a global property, since in order to check it one has to consider 
the relation between the given set of vertices and the entire graph. In this sense, 
determining the domatic number of a graph intuitively appears to be harder than 
computing its chromatic number, notwithstanding that both problems are NP-complete. 
More to the point, the algorithms developed for graph coloring seem to be of no help 
in designing algorithms for dominating set or domatic number problems. 

After introducing some definitions and notation in Section [2] we describe and 
analyze our algorithm in Section [3J the actual pseudo-code is shifted to the appendix. 
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In Section |4] we give another deterministic and a randomized algorithm, which have 
an even better running time for graphs with small maximum degree. Finally, we 
summarize and discuss our results in Section|5] 

2 Preliminaries and Simple Observations 

We start by introducing some graph-theoretical notation. We only consider simple, 
undirected graphs without loops in this paper. Let G = (V, E) be a graph. Unless stated 
otherwise, n denotes the number of vertices in G. The neighborhood of a vertex v in 
V is defined by N(v) = {u E V | {u, v} E E}, and the closed neighborhood of v 
is defined by N[v] — N(v) U {v}. For any subset S C V of the vertices of G, 
define N[S] = \J veS N[v] and N(S) = N[S] - S. The degree of a vertex v in G 
is the number of vertices adjacent to v, i.e., deg G (v) — \\N(v)\\. If the graph G is 
clear from the context, we omit the subscript G. Define the minimum degree in G 
by min-deg(G) = mm v£ v deg(v), and the maximum degree in G by max-deg(G) = 
maxugv' deg(v). A path Pk = u\Ui ■ ■ ■ Uk of length k is a sequence of k vertices, 
where each vertex is adjacent to its successor, i.e., {m, u^+i} E E for 1 < i < k — 1. 
If, in addition, {iik, ui} £ E, then path Pk is said to be a cycle, and we write C% 
instead of Pk. 

Definition 1 Let G — (V, E) be a graph. A subset D C V is a dominating set of G if 
and only if N[D] = V, i.e., if and only if every vertex in G either belongs to D or has 
some neighbor in D. The domination number of G, denoted 7(G), is the minimum size 
of a dominating set of G. The domatic number of G, denoted 5(G), is the maximum 
number of disjoint dominating sets of G, i.e., (5(G) is the maximum k such that V — 
V\ U V% U . . . U Vk, where ViHVj = §for 1 < i < j < k, and each Vi is a dominating 
set of G. The dominating set problem asks, given a graph G and a positive integer k, 
whether or not 7(G) < k. The domatic number problem asks, given a graph G and a 
positive integer k, whether or not 5(G) > k. 

For fixed k > 3, both the dominating set problem and the domatic number problem 
are known to be NP-complete, see Garey and Johnson IGJ79I . Thus, they are not 
solvable in deterministic polynomial time unless P = NP, and all we can hope 
for is to design an exponential-time algorithm having a better running time than the 
trivial exponential time bound. For exponential-time algorithms, it is common to drop 
polynomial factors, as indicated by the O notation: For functions / and g, we write 
f E O(g) if and only if / E 0(p ■ g) for some polynomial p. The naive deterministic 
algorithm for the dominating set problem runs in time 0(2 n ). Fomin, Kratsch, and 
Woeginger [FKW04| improved this trivial upper bound to 0(1.93782"). For various 
restricted graph classes, they achieve even better bounds. 

The naive deterministic algorithm for the domatic number problem works 
as follows: Given a graph G and an integer k, it sequentially checks every 
potential solution (i.e., every possible partition of the vertex set of G into k sets 
Di, D2, ■ ■ ■ , Dk), and accepts if and only if a correct solution is found (i.e., if and 
only if each Di is a dominating set). How many potential solutions are there? The 
number of ways of partitioning a set with n elements into k nonempty, disjoint 
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subsets can be calculated by the Stirling number of the second kind: S^n, k) — 

H Ej=o 1 ( _1 ) i (j)( fc ~ *)"' which y ields a runnm g time of 6(k n ). A better result 
can be achieved via the dynamic programming across the subsets technique, which 
was introduced by Lawler |Law76| to compute the chromatic number of a graph 
by exploiting the fact that every minimum chromatic partition contains at least one 
maximum independent set. By suitably modifying this technique, one can compute the 
domatic number of a graph in time (9(3™). This is done by generating all dominating 
sets of the graph with increasing cardinality, which takes time 

(?V = (1 + 2)™ = 3". 

The difference to Lawler' s algorithm lies in the fact that all dominating sets need 
to be checked, whereas only maximum independent sets are relevant to compute the 
chromatic number. 

Proposition 2 Let G = (V, E) be a graph. Then, the domatic number 5(G) can be 
computed in time (9(3"). 

One tempting way of designing an improved algorithm for the domatic number 
problem might be to exploit the result for the dominating set problem mentioned above. 
However, we observe that no such useful connection between the two problems exists 
in general. The first part of Proposition [5] shows that an arbitrary given minimum 
dominating set is not necessarily part of a partition into a maximum number of 
dominating sets. The second part of Proposition [5] shows that, given an arbitrary 
partition into a maximum number of dominating sets, it is not necessarily the case 
that one set of the partition indeed is a minimum dominating set. Thus, for solving the 
domatic number problem, one cannot use in any obvious way the exact (5(1.93782™) 
algorithm for the dominating set problem by Fomin et al. |FKW04|. Proposition|5]is 
stated for graphs with domatic number 3; it can easily be generalized to graphs with 
domatic number k > 3. The proof of Proposition[5]can be found in the appendix. 

Proposition 3 1. There exists some graph G with 5(G) — 3 such that some 
minimum dominating set D ofG is not part of any partition into three dominating 
sets of G. 

2. There exists some graph H = (V 1 E) with 5(H) = 3 such that for each partition 
V = D1UD2UD3 into three dominating sets of H and for each i, \ \Di\ \ > "/(H). 

For the three domatic number problem, no algorithm with a running time better 
than (9(3™) is known. We improve this trivial upper bound to (9(2.9416"). 

We now define some technical notions suitable to measure how "useful" a vertex is 
to achieve domination of the graph G — (V, E). Intuitively, the vertex degree is a good 
(local) measure, since the larger the neighborhood of a vertex is, the more vertices are 
potentially dominated by the set to which it belongs. The technical notions introduced 
in Definition 0] will be used later on to describe our algorithm. 
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Definition 4 Let G = (V, E) be a graph with n vertices, and let V = (Di, D2, D%, R) 
be a partition of V into four sets, D\, D2, D3, and R. The subsets Di of V will 
eventually yield a partition of V into the three dominating sets (if they exist) to be 
constructed, and the subset R C V collects the remaining vertices not yet assigned at 
the current point in the computation of the algorithm. Let r = \\R\\be the number of 
these remaining vertices, and let d = n — r be the number of vertices already assigned 
to some set Di. The area of G covered by V is defined as area-p (G) = Y^=i I l-M-^i] 1 1- 
Note that area-p(G) = 3n if and only if D\, D2, and D3 are dominating sets of G. For 
a partition V, we also define the surplus of graph G as surplus.p(G) = area^(G) — 3d. 

Some of the vertices in R may be assigned to three, not necessarily disjoint, 
auxiliary sets A\, A2, and A3 arbitrarily. Let A = (A\, A2, A3). For each vertex 
v e R and for each i with 1 < i < 3, define the gap of vertex v with respect to set Di 
by 

&VvA v ^)-\ j_ otherwise, 

where _L is a special symbol that indicates that gap^, . (v, i) is undefined for this v 
and i. ( Our algorithm will make sure to properly handle the cases of undefined gaps.) 
Additionally, given V and A, define for all vertices v G _R: 

maxgapp A (v) = max{gap 7 , <A (v, i) \ 1 < i < 3}, 
mingapp A (v) = mm{gap VA (v,i)\l<i<3}, 

3 

sumg&v VA (v) = ^2gap VA {v,i). 

i=l 

Given G, V, and A, define the maximum gap of G and the minimum gap of G by 
taking the maximum and minimum gaps over all vertices in G not yet assigned: 

maxgapp A (G) = max{maxgap p A (v) \ v G R}, 
mingapp <A (G) = min{mingap-p <A (v) | v G R}. 

Let V be given. A vertex u G V is called an open neighbor of v G V if u G N[v] 
and u has not been assigned to any set D\, D 2 , or D 3 yet. A potential dominating 
set Di, 1 < i < 3, is called an open set of v G V if its closed neighborhood does 
not include v, i.e., v is not dominated by D^ The balance of v G V is defined as the 
difference between the number of open vertices and the number of open sets. Formally, 
define 

openNeighborSp (v) = {u G N[v] \ u G R}, 

openSetSp(u) = {i e {1, 2, 3} | v (£ iV[A]}, 
balancep(u) = | |openNeighborSp(u)| | — | |openSets T ,(u)| |. 

We call a vertex v G V critical if and only if balancep(w) < and ||openSets T ,(u)|| > 
0. 
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The proof of the next proposition is straightforward. Once balance-p(v) = 0, no 
two vertices remaining in N[v] n R can be assigned to the same dominating set Di, 
1 < i < 3, since balance-p(w) would then be negative. 

Proposition 5 Let V = {D\, D2, D3, R) be given as in Definition^, and v £ V be 
a critical vertex for this partition. The only way to modify V so as to contain three 
dominating sets is to assign all vertices u £ N[v] D Rto distinct dominating sets Di. 

3 The Algorithm 

Our strategy is to recursively assign the vertices v £ V to obtain a correct potential 
solution consisting of a partition into three dominating sets, D\, D2, and D3. Once a 
previous assignment of v to some set Di turns out to be wrong, we remember this by 
adding this vertex to Aj. More precisely, the basic idea is to first pick those vertices 
with the highest maximum gap. While the algorithm is progressing, it dynamically 
updates the gaps for every vertex in each step. We now state our main result. 

Theorem 6 The three domatic number problem can be solved by a deterministic 
algorithm running in time 0(2.9416"). 

Proof. Let G = (V, E) be the given graph. The algorithm seeks to find a partition 
of V into three disjoint dominating sets. Note that every vertex v £ V is contained in 
one of these sets and is dominated by the remaining two sets, i.e., it is adjacent to at 
least one of their elements. The algorithm is described in pseudo-code in the appendix, 
see Figures 13 13] EE] |5] and[6] Since 6(G) < min-deg(G) + 1, we may assume that 
min-deg(G) > 2. 

The algorithm starts by initializing the potential dominating sets D±, D2, and D3 
and the auxiliary sets A%, A2, and A3, setting each to the empty set. The initial partition 
thus is V = (0, 0, 0, V) and the initial triple of auxiliary sets is A = (0,0,0). 

Then, the recursive function DOMINATE is called for the first time. It is always 
invoked with graph G, a partition V = (Di, D2, D3, R), and a triple A = (^1,^2,^3) 
of not necessarily disjoint auxiliary sets. V and A represent a situation in which the 
vertices in V — R have been assigned to D\, D2, and D3, and v £ A{ means that 
in some previous recursive call to function DOMINATE the vertex v has been assigned 
to Di without successfully changing V to contain three dominating sets. 

Function DOMINATE starts by calling Recalculate-Gaps, which calculates all 
gaps with respect to V and A. Additionally, openNeighbors T ,(u), openSetSp(z;), and 
balance^(w) are determined for every vertex v £ V. Four trivial cases can occur. 

Case 1: The sets D\, D2, and D3 are dominating sets of graph G. In this case, we are 
done and may add the remaining vertices v £ R to any set Di, say to D\. 

Case 2: For some vertex v £ V, we have balance^ (v) < 0. That is, there are less 
vertices in R n N[v] than dominating sets with v ^ N[Dj\. Thus, no matter how the 
vertices in R n N[v] are assigned, V won't contain three dominating sets. We have 
run into a dead-end and return to the previous level of the recursion. 
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Case 3: There exists a vertex v £ R that is also a member of two of the auxiliary sets 
A\, A2, and A 3 . Hence, vertex v was previously assigned to two distinct sets Di 
and Dj, 1 < i < j < 3, but the recursion returned without success. We assign v to 
the only possible set Dk left, with i ^ k ^ j. 

Case 4: For some vertex v £ V", we have balance-p(w) = and | |openSetSp(u)| | > 0. 
That is, v is a critical vertex, since it is not dominated by all three sets D\, D 2 , and 
D 3 contained in the current V, and there are as many open neighbors as open sets left 
for it. Note that this is the case for each vertex v with deg(v) — 2 and N[v] n R 7^ 0, 
as v and its two neighbors have to be assigned to three different dominating sets. 
We select one of the at most three vertices left in N[v] n R, say u, and call function 
Assign(G,V,A, u, i) for all i with u <£ A { . 

Function Handle-Critical- Vertex deals with the latter three of these trivial 
cases. After they have been ruled out, one of the remaining vertices v £ R is selected 
and assigned to one of the three sets Dj, under the constraint that a vertex v £ R cannot 
be added to Di if it is already a member of A t . This case occurs whenever the recursion 
returns because no three dominating sets could be found with this combination. The 
recursion continues by calling Assign(G, V, A, v, i), which adds v to D i7 and then 
calls Dominate(G,7 : ',^4). If no three dominating sets are found by this choice, we 
remember this by adding v to the set At. A final call to DOMINATE is made without 
assigning a vertex to one potential dominating set Di. If this call fails, the recursion 
returns to the previous level. This completes the description of the algorithm. We now 
argue that it is correct and estimate its running time. 

To see that the algorithm works correctly, note that it outputs three sets D\, D 2 , 
and D 3 only if they each are dominating sets of G. It remains to prove that these sets 
are definitely found in the recursion tree. All drop-backs within the recursion occur 
when, for the current V = (D\, D2, -D3, R), we have balance-p(u) < for some 
vertex v e V. Thus, V cannot be modified so as to contain a correct partition into 
three dominating sets on this branch of the recursion tree. Since the algorithm checks 
every possible partition of G into three sets, unless it is stopped by such a drop-back, 
a partition into three dominating sets will be found, if it exists. If the algorithm does 
not find three dominating sets, it eventually terminates when returning from the first 
recursive call of function DOMINATE. It reports the failure, and thus always yields the 
correct output. 

To estimate the running time of the algorithm, an important observation is that the 
recalculation of the gaps takes no more than quadratic time in n, the number of vertices 
of the graph G. Thus, in terms of the 0-notafion, the running time of the algorithm 
depends solely on the number of recursive calls. Let T(m) be the number of steps of 
the algorithm, where m is the number of potential dominating sets left for all vertices 
that have not been selected as yet. Initially, every vertex may be a member of any of 
the three dominating sets to be constructed (if they exist), hence m = 3n. 

There are two scenarios where the algorithm calls function DOMINATE recursively. 
If Handle-Critical-Vertex detects a vertex v <E V as being critical, it selects a 
vertex u e N[v] n R and calls function Assign (and thus Dominate) for each i with 
u Ai. Since every critical vertex v £ V remains critical as long as N[v] n R ^ 0, 
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function Handle-Critical-Vertex will be called until all vertices in N[v] n R 
have been assigned to any of D\, D%, and D3. Since ||openSets-p(ii)|| < 3, at most 
three vertices in the closed neighborhood of v have not been assigned when v turns 
critical. By Proposition [5] all vertices in N[v] D R have to be assigned to different 
dominating sets. If ||openNeighborSp(w)|| = 3, we have at most six combinations; if 
we have two open neighbors for a critical vertex, there are at most two combinations 
left; and finally, for one open neighbor u € N[v] n R, there remains only one possible 
choice to assign u to one of the sets D\, D2, and D3. Thus, in the worst case, we have 
T{m) < QT(m — 6), as we will handle three vertices for which at least two choices 
for dominating sets are left. With m = 3ri, it follows that T{m) < 6 m / 6 = 6"/ 2 , i.e., 
T(m) = 0(2.4495"). 

The only other branching into two different recursive calls happens in the main 
body of function DOMINATE, when selecting a vertex v with the currently highest 
maximum gap with respect to V and A. Two cases might occur. On the one hand, we 
might have considered a correct dominating set Di for v. If v had not been looked at 
so far, i.e., if v is not contained in any set Aj, 1 < j < 3, j 7^ i, we have eliminated 
all three possible sets for v to belong to. Thus, in this case, T(to) — T(m — 3). On 
the other hand, if the algorithm returns from the recursion and thus did not make the 
right choice for v, we have T{m) = T(m — 1), since v is added to A i7 and function 
Dominate is called without assigning any vertex. Summing up, we have T(m) < 
T(m — 1) + T(m — 3). In the second case, we have already visited vertex v in a 
previous stage of the algorithm and unsuccessfully tried to assign it to some set Dj, 
with 1 < j < 3. There are only two dominating sets for v left. Either way, if we 
put v into the correct dominating set right away or fail the first time, we have T(m) = 
T(m — 2). Summing up both cases, we have T(m) < 2T(m — 2). Suppose that the 
first and the second case occur equally often, i.e., the algorithm considers every vertex 
twice. It then follows that 

T(m) < -(T(m- 1) +T(m-3)) + -(2T(m- 2)) 

with to = 3rt. Thus, we have T(m) = (D(3 n ), and the trivial time bound cannot be 
beaten. To improve this running time, we have to make sure that the recursion tree 
will not reach its full depth, i.e., not all vertices are considered by the algorithm or 
function Handle-Critical-Vertex will be called for a sufficiently large portion 
of the vertices. It is clear that the algorithm has found three dominating sets once 
areap(G) = in (recall the notions from Definition 0J. By selecting the maximum 
gap possible for a partition V, we try to reach this goal as fast as possible. For every 
vertex v E R that we assign to one of the potential dominating sets Di, 1 < i < 3, 
we increase areap(G) by gap p A (v, i), and additionally we add (gap v A (v, i) — 3) to 
surplus T ,(G). 

Since the vertices of degree two are critical, they and their neighbors can be handled 
in time O (2.4495™), as argued above. So assume that min-deg{G) > 3. Then, we have 
maxgapp A (G) > 3 at the start of the algorithm. If this condition remains to hold for 
at least 3rt/4 steps, we have reached area-p(G) = in, and the algorithm terminates 
successfully. To make use of more than 3n/4 vertices, maxgap p A (G) has to drop 
below four at one point of the computation. We exploit the fact that up to this point, 
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the surplus has grown sufficiently large with respect to n. Decreasing it will force 
maxgapp A (G) to drop below three, and this condition can hold only for a certain 
portion of the remaining vertices until the algorithm terminates. To see this, we now 
analyze the remaining steps of the algorithm after the given graph G has reached a 
certain maximum gap with respect to the current V and A. 

If maxgapp A (G) = 0, the recursion stops immediately. Either we have already 
found three disjoint dominating sets (in which case we put the remaining vertices v E R 
into set D\ and halt), or one vertex has not been dominated by one set Di in V yet. 
Since no positive gaps exist for the vertices v S R, V cannot be modified to a valid 
partition into three dominating sets. Function Handle-Critical-Vertex returns 
true immediately after detecting balancep(i;) < for some vertex v G V, and function 
Dominate drops back one recursion level. The question is how many vertices are left 
in R when we reach maxgapp A {G) = 0. 

Lemma 7 Let G = (V, E) be a graph and V = [D\, -D3, R) be a partition ofV 
as in Definition^ Let r = \ \R\ \ and maxgapp A (G) = 3. Then, for at least r / '64 
vertices in R, the algorithm will not recursively call function DOMINATE. 

Proof of Lemma El Let maxgapp ^(G) = k with k > 0. Since gap VA (v, i) < k 
for each v E R and for each i, 1 < i < 3, we have J2 v eR sum g a P-p a( v ) — ^kr. 
Every vertex v that is selected for a set Di with gap.p A (v, i) — k decreases at least k 
gaps of the vertices in R — {v} by one. Otherwise, Handle-Critical-Vertex 
would have found a critical vertex u E N[v] with N[u] n R — {v}. Then, 
either ||openSets.p(w)|| > 1 (which implies balancep(u) < and we abort), or 
I |openSets-p(u) 1 1 = 1, in which case v is added to the appropriate set Di without further 
branching of function Dominate. Thus, if no critical vertex is detected, selecting a 
vertex v E R for some set Di decreases at least k gaps, and since v does not belong 
to R anymore, additionally all gaps previously defined for v are now undefined. So 
the lowest possible rate at which the gaps are decreased is related to the maximum gap 
ofG 

Now suppose that maxgap^, A {G) = 3 and sumgap^, ,(v) = 9 for all vertices 
v E R. We always select a vertex v with the highest summation gap of all vertices 
u E R with maxgap-p ^(u) = 3. As long as there exists a vertex v E R with 
gapp A (v, i) = 3 for all i, it will be selected by the algorithm. After calling function 
RECALCULATE-GAPS, the number of gaps equal to three will be decreased at least 
by six. If exactly three other gaps of vertices in R — {v} decrease by one in every 
step, it takes at least r/4 vertices until sumgapp A (v) < 9 for all v E R. Another 
1/4 of the 3r/4 vertices remaining have to be selected until sumgapp_ 4 (f) < 8. 
Adding 1/4 of the 9r/16 vertices left in R, we have reached maxgapp A (G) = 2 
with sumgapp ,(v) = 6 for all vertices v E R. This implies that every defined gap is 
equal to two. Summing up, we have selected 

1 13 19 37 
— • r H — • — r H — • — r = — r 

4 4 4 4 16 64 

vertices until maxgapp A (G) = 2, under the constraint that a minimum number of gaps 
is reduced in each step, while simultaneously trying to reduce the maximum summation 
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gap in the fastest possible way. This way we reach level maxgap^, A (G) = with as 
few vertices left in R as possible, which describes the worst case that might happen. 

Analogously, we can show that maxgap p A (G) drops from 2 to 1 after selecting 
another 19r/64 vertices. And once we have maxgap p A (G) = 1, it takes 7r/64 
vertices to get to maxgap p A (G) = 0. Now, there are r/64 vertices remaining in R, 
which do not have to be processed recursively. | Lemma[7] 

Continuing the proof of Theorem |6] note that we assumed min-deg(G) > 3, so 
when the gaps are initialized for graph G, we have mingap p A (v) > 4 for each vertex 
v € V. Thus, more than three vertices are dominated by the selected set Dj for 
vertex v. As long as maxgap-p A {G) > 3 is true, surplus-p(G) is increasing. The only 
way to lower the surplus is by adding vertices v to a set Di with gap-p A (v, i) < 3. 
The surplus decreases by one when gap-p A (v, i) = 2, and it decreases by two when 
gap VA {v,i) = 1. 

Let S = surplus,, (G) be the surplus collected for a partition V until we reach a 
point where maxgapp A (G) = 3. To make use of the most recursive calls and to even 
out the surplus completely, there have to be at least r — \ \R\ \ vertices remaining with 

„ 37r „ 19r „ 7r 
°-lM +1 -lM +2 -64^ 

so r > 64S733. A fraction of 1/64 of these vertices will be handled by the algorithm 
without branching into more than one recursive call, which is at least 5/33. The 
question is how big the surplus S might grow and how many vertices are left in R 
before maxgapp^G) = 3 is reached. The lowest surplus with as few vertices in 
R as possible occurs if min-deg(G) = max-deg(G) = 3. Surplus S is increased by 
one in each step until we arrive at maxgap^, A (G) = 3. When selecting a vertex 
v of degree 3 for a set D i7 the gap of its neighbors u G N(v) and the gaps of the 
neighbors of every u might be decreased. Summing up, at most 1 + 3 + 3- 2 = 10 
vertices can have decreased their gaps for some i. After selecting at least n/10 vertices 
for each i, we have mingap^, A (G) — 3 (in the worst case). From this point on, we 
cannot be sure if the next vertex selected for some Di satisfies gap-p A (v,i) > 3. 
But so far we have already collected a surplus of S = 3n/10, and applying this we 
obtain 64n/110 < r < 7n/10. Thus, for at least n/110 vertices we never branch into 
two different recursive calls. Setting m = 3(109n/110), we obtain a running time 
of 0(2.9416"). | 



4 Graphs with Bounded Maximum Degree 

As seen in the last section, the running time of the algorithm crucially depends on 
the degrees of the vertices of G. If we restrict ourselves to graphs G with bounded 
maximum degree (say A = max-deg(G)), we can optimize our strategy in finding three 
disjoint dominating sets. In this section, we present a simple deterministic algorithm, 
which has a better running time than the algorithm from Theorem[6] provided that A 
is low. By using randomization, we can further improve the running time for graphs G 
with low maximum degree. 
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Before stating the two results, note that graphs with maximum degree two can 
trivially be partitioned into three dominating sets, if such a partition exists. Every 
component of such a graph is either an isolated vertex, a path, or a cycle, and each such 
property can be recognized in polynomial time. 

Proposition 8 Let G — (V, E) be a given graph with max-deg(G) = 2. There exists a 
partition of the vertices of G into three dominating sets if and only if every component 
of G is a cycle of length k such that 3 divides k. 

We use the terms from Definition 0] in Section[3]to describe a snapshot within the 
algorithm. For any partition V = (Dx, D2, -D3, R), some vertices of V have already 
been assigned to the potential dominating sets Dx, D2, and D3, while all the remaining 
vertices are in R. The auxiliary sets A = (A\, A2, A3) will not be needed in this 
section. Only connected graphs are considered, as it is possible to treat every connected 
component separately, producing the desired output within the same time bounds. 

TableQJlists the running times of both the deterministic and the random algorithm, 
where the maximum degree of the input graph is bounded by A, 3 < A < 8. Note that 
the exact deterministic algorithm from Theorem|5]in Section[3]beats the deterministic 
algorithm from Theorem|9]whenever A > 7. 



A 


3 


4 


5 


6 


7 


8 


deterministic 


2.2894" 


2.6591" 


2.8252" 


2.9058" 


2.9473" 


2.9697" 


randomized 


2" 


2.3570" 


2.5820" 


2.7262" 


2.8197" 


2.8808" 



Table 1: Results for max-deg(G) = k, where 3 < k < 8 



Theorem 9 Let G = (V, E) be a graph with max-deg(G) — A, where A > 3. 
There exists a deterministic algorithm solving the three domatic number problem in 
time 0(d~z), where 

A-2 



A\ /A 



E 

6=1 



(4.1) 



Proof. The algorithm works as follows. We start with an arbitrary vertex v 6 V and 
assign it to the first set D\. In each step, we first check whether we found a partition 
V = (D±, D2, D3, R) into dominating sets D\, D2, and D3. If not, one vertex v € V 
is selected that is not dominated by all three sets D\, D2, and D3, and additionally 
has a vertex u E N[v] in its closed neighborhood that has already been added to some 
set D h 1 < i < 3. It follows that 1 < ||openSets P (i;)|| < 2. 

If balancep(w) < 0, we return within the recursion. Otherwise, we try all 
combinations to partition the vertices in N[v] n R, so that after this step vertex v is 
dominated by all three potential dominating sets. If no such combination leads to a 
valid partition, we again return within the recursion. 

Suppose now that balance-p (v) > 0, ||openSetSp(w)|| = 2, and N[v] flfli ^ 0. 
To obtain three disjoint dominating sets, at least one vertex in N[v] has to be assigned 
to D2, and at least one vertex in N[v] has to be added to D3. This limits our choices, 
especially if the degree of v is bounded by some constant A. 
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To measure the running time of the algorithm, we consider the worst case with 
the most possible combinations that might yield a partition into three dominating sets. 
This occurs when only one vertex u £ N[v] has already been added to one set, i.e., 
\\N[v] n (Di UD 2 U £> 3 )|| = 1. lfN[v] n Dx ^ 0, then any number between and 
A — 2 of vertices in N [v] n R may be assigned to the same set D\ . Let this number be a. 
It follows that from one to A — a — 1 vertices remaining in N[v] n R are allowed to be in 
the next potential dominating set This is how Eauation l4.1l for d is derived. After 
assigning the last vertices in N[v]DR to the dominating set D3, exactly A vertices have 
been removed from R. Thus, we have a worst case running time of 0{d~s). Tabled 
lists the running time for graphs with maximum degree from three to nine. I 

In the next theorem, randomization is used to speed up this procedure. Instead of 
assigning all vertices in the closed neighborhood of some vertex v G V in one step, 
only one or two vertices in N[v] OR are added to the potential dominating sets D±, 
D2, and D3. The goal is to dominate one vertex by all three sets in one step. We will 
select the one or two vertices that are missing for this goal at random. 

Theorem 10 Let G — (V, E) be a graph with max-deg(G) — A, where A > 3, and 
let d be defined as in Equation i4.lt in Theorem |P| For each constant c > 0, there 
exists a randomized algorithm solving the three domatic number problem with error 
probability at most e~ c in time 0(r*), where 

A 

(4.2) 



3 A-2- 

Proof. Let graph G = (V, E) be given with max-deg(G) = k. As in the 
deterministic algorithm, we start by adding a random vertex to the set D\. In every 
following step, a vertex v £ V is selected with < | |openSets T ,(u)| | < 3, so it is 
N[v] n (D 1 UD 2 U D 3 ) ^ 0. If ||openSets.p(u)|| = 1, we have N[v] n A = for 
one i with 1 < i < 3. We randomly choose a vertex we N[v] n R and assign it to set 
Di, in order that v is dominated by all three sets afterwards. If | |openSetSp(u)| | = 2, 
we randomly select two vertices u\, 112 € R in the closed neighborhood of v. Another 
random choice is made when deciding how to distribute these two vertices among the 
two potential dominating sets that have not dominated v up to now. 

Suppose G indeed has a partition into three dominating sets. We have to measure 
the error rate when making our random choices to estimate the success probability 
of the algorithm. In every step, a vertex v £ V is selected with at least one vertex 
u G N[v] in its closed neighborhood that has already been added to one of the sets D\, 
D2, or D3 . The highest error occurs when exactly one vertex in N[v] is not included 
in R, so we restrict our analysis to this case. To obtain a valid partition into three 
dominating sets, there are at most d choices left to partition the vertices remaining in 
N[v] n R. Here, d is the number from Eauation l4.il Once we selected and assigned 
two vertices from N[v] n R at random, there are 3 fc_2 possibilities left to partition the 
vertices in the closed neighborhood of v that are still left in R. Our success rate when 
selecting the two vertices is therefore 3 fe ~ 2 /d. 

To achieve an error probability of below e~ c , the algorithm needs to be executed 
more than once. The repetition number of the algorithm equals the reciprocal of the 
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success rate, which explains Equation 14.21 Since two vertices are processed in every 
step, the overall running time is 0(r% ). I 



5 Conclusion 

We have shown that the three domatic number problem can be solved by a deterministic 
algorithm in time (5(2.9416"). Furthermore, we presented two algorithms solving the 
three domatic number problem for graphs with bounded maximum degree, improving 
the above time bound for graphs with small maximum degree. Although our running 
times seem to be not too big of an improvement of the trivial 0(3" ) bound, they 
are to our knowledge the first such algorithms breaking this barrier. For k > 3, the 
decision problem of whether 8(G) > k can be solved in time 0(3") by Lawler's 
dynamic programming algorithm for the chromatic number, appropriately modified for 
the domatic number problem. Therefore, it would not be reasonable to use our gap 
approach of Section[3]to decide if 8(G) > k for a graph G and k > 3. 

Acknowledgement. We thank Dieter Kratsch for pointing us to Lawler's algorithm. 
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A Proof of Proposition |3] 



Proof. Figure ^ shows the graphs G and H whose existence is claimed. In this 
figure, the numbers i\j within a vertex have the following meaning: i indicates which 
dominating set Di this vertex belongs to in a fixed partition into three dominating sets, 
and j indicates a specific choice of a minimum dominating set S of the graph by setting 
j = 1 if and only if this vertex belongs to S. 




Figure 1 : Graphs G and H for Proposition|5] 

For the first assertion, look at the graph G shown on the left-hand side of Figure^ 
Note that 7(G) = 2. In particular, D = {u 3 , 1/5} is a minimum dominating set of G. 
Note further that 5(G) — 3. In particular, a partition into three dominating sets of 
G is given by D\ = {m.u^uj}, D 2 = {u2,u 5 }, and D 3 = {u 3 ,uq}. However, 
D cannot be part of any partition into three dominating sets, since the only neighbors 
of U4, namely 113 and M5, belong to D. 

Note that the minimum dominating set D 2 = {1*2, u s} of G defined above indeed 
is part of a partition into three dominating sets. The second part of the proposition, 
however, shows that this is not always the case. Consider the graph H = (V, E) shown 
on the right-hand side of Figure ^ We have j(H) = 2 by choosing the minimum 
dominating set D = {1)1,1)2}, which is unique in this case. Again, 8(H) = 3. The 
only way, up to isomorphism, to partition the vertex set of H into three dominating 
sets is given by D x = {vi,v 7 ,v&}, D 2 = {v 2 , v 6 , v 9 }, and D 3 = {v 3 ,v 4 ,v 5 }. Thus, 
min{||Di||, H-D2II, H-D3II} > l(H) for each partition into three dominating sets. I 
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B Pseudo-Code of the Algorithm from Theorem |6] 

Figures[2][3]|4]|5] and[6]describe the algorithm from Theorem[6]in pseudo-code. 



Algorithm for the Three Domatic Number Problem 

Input: Graph G — (V, E) with vertex set V = {v±, V2, ■ ■ ■ , v n } and edge set E 
Output: Partition of V into three dominating sets D\, D2, -D3 C V or "failure" 

Set each of D\, D2, D3, A\, A2, and A3 to the empty set; 
Set R = V; 

SetV =(D 1 ,D a ,D 3 ,Ry, 
SetA=(Ai,A 2 ,A 3 ); 

Dominate(G, V, A); // Start recursion 

output "failure" and halt; 

Figure 2: Algorithm for the Three Domatic Number Problem 
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Function Dominate(G, V, A) { // V is a partition of graph G, 

II A is a triple of auxiliary sets 

Recalculate-Gaps(G, V, A); 
if (each Di is a dominating set) { 

£>i = £>i U R; 

output D 1 ,D 2 ,D 3 ; 

} 

if ( not Handle-Critical- Vertex(G, V, A)) { 
select vertex v € R with 

maxgapp A (v) = maxgap^, A (G) and 

sumgapp A (v) = maxjsumgapp A (u) \ u e R A rnaxgap^, ^(m) = 
maxgapp ^(G)}; 
find « with gap.p A (v, i) = maxgap,p A (v); 

Assign(G, V, A, v, i); //First recursive call 

Ai = AiU {v}; 1 1 If recursion fails, put v in Ai and try again 

D0MINATE(G, V, A); // Second recursive call 

} 

return; 

} 



Figure 3: Recursive function to dominate graph G 



Function Assign(G, V, A, v, i) { 

Di = Di U {v}; 
R = R-{v}; 
DoMiNATE(G,7 : ',yl); 

} 

Figure 4: Function to assign vertex v to set Di 



18 



Function Recalculate-Gaps(G, V, A) { //Visa partition of graph G, 

1/ A is a triple of auxiliary sets 

for all (vertices v € V) { 
if (vertex v e R) { 
for all (i = 1,2,3) { 

if (v £ A t ){gw VA {v,i) = \\N[v]\\ - \\{u G JV[w] | (3«; G G 

A]}||;} 

else gapp ^(u, i) — _L ; // _L indicates that gap P ^(w, i) is undefined 

} 

maxgapp iiA (u) = max ie{ i i2 ,3}{gap Pi ^(u, i)}; 
mingapp iiA (u) = min 4(E{1 , 2 , 3 }{gap P >iA (u, «)}; 
sumgapp = £\ £{1 2 3} gap P)iA (u, i); 

} 

openNeighbors-p(v) ={ii£ 7V[v] | u G i?}; 
openSetSp(u) = {ie {1, 2, 3} | v £ N[D t ]}; 
balance-p(w) = ||openNeighbors 7 ,(w)|| — ||openSetSp(t;)||; 

} 

maxgapp iA (G) = max veR {maxgap-p A (v)}; 
mingapp A (G) = min„ eiJ {rningap.p tA (v)}; 



Figure 5: Function to recalculate gaps after partition has changed 



19 



Function boolean Handle-Critical- Vertex(G, V, A) { 

for all (vertices v e V) { 

if (balance-p(w) < 0) { // impossible to three dominate v with V 

return true; 

} else if (||{i e {1,2, 3} | v e Ai}\\ == 2) { // one choice for v remaining 
select i with v Af, 
ASSlGN(G,V,A,v,i); 
return true; 

} else if (balance-p(i;) == and ||openSets P (t;)|| > 0) { //vis critical 

select u e N[v]C\ R; 

for all (i with u £ Ai and v not dominated by DA 

Assign(G,7 : ',^,u, i)\ 
return true; 

} 

} 

return false; // no critical vertices were found 

} 

Figure 6: Function to handle the critical vertices 
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